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ABSTRACT 

Context. Globular clusters are an excellent laboratory for stellar population and dynamical research. Recent studies have shown that 
these stellar systems are not as simple as previously assumed. With multiple stellar populations as well as outer rotation and mass 
segregation they turn out to exhibit high complexity. This includes intermediate-mass black holes which are proposed to sit at the 
centers of some massive globular clusters. Today's high angular resolution ground based spectrographs allow velocity-dispersion 
measurements at a spatial resolution comparable to the radius of influence for plausible IMBH masses, and to detect changes in the 
inner velocity-dispersion profile. Together with high quality photometric data from HST, it is possible to constrain black-hole masses 
by their kinematic signatures. 

Aims. We determine the central velocity-dispersion profile of the globular cluster NGC 2808 using VLT/FLAMES spectroscopy. In 
combination with HST/ ACS data our goal is to probe whether this massive cluster hosts an intermediate-mass black hole at its center 
and constrain the cluster mass to light ratio as well as its total mass. 

Methods. We derive a velocity-dispersion profile from integral field spectroscopy in the center and Fabry Perot data for larger radii. 
High resolution HST data are used to obtain the surface brightness profile. Together, these data sets are compared to dynamical models 
with varying parameters such as mass to light ratio profiles and black-hole masses. 

Results. Using analytical Jeans models in combination with variable M/L profiles from N-body simulations we find that the best 
fit model is a no black hole solution. After applying various Monte Carlo simulations to estimate the uncertainties, we derive an 
upper limit of the back hole mass of Mbh < 1 x 10"* Mq (with 95% confidence limits) and a global mass-to-light ratio of M/Ly = 
(2.1 +0.2) Mo/Lq. 

Key words, black hole physics - globular cluster: individual (NGC 2808) - stars: kinematics and dynamics 
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1. Introduction young star clusters by runaway merging. However, taking 

into account mass loss through stellar winds, lYungelson et al.l 

^ ; Kinematics of globular clusters have long been a field of interest ,2008) found that super-massive stars with initial masses up 
L; . in observational and computational astronomy. In the last years, ^qOO Mq reduce to objects less massive than ~ 150 Mq 

. ^ ■ a new aspect, seai'chmg for the signature of central mtermediate- ^y the end of their l ives. An other scenario was presented 

^ mass black holes (IMBHs) was added to the picture. The spe- ^y Miller & Hamilto?] dSh who discussed the formation of 

H cific mass range of black holes (10 - 10 Mo) had not been ob- intermediate-mass black holes through mergers of stellar-mass 

. ? ; served before and for a long time intermediate-mass black holes ^lack holes in globular clusters when starting with a M > 50 Mo 

even thought not to exist. However, recent observations (eg. gggd ^i^ck hole. Further, they presented scenarios for the capture 



Pebhardt et al. 2000; Gerssen et al. 2002; Gebhardt et al. 2005; dusters by their host galaxies and accretion in the galactic 
i Noyola et al. 2008; Lutzgendorf et al. 201 1) have shown that the ^jj^k in order to explain the obsei-ved bright X-ray sources. In 



velocity-dispersion profiles of some globular clusters and dwarf addition, remnants of massive Population III stars could have 
galaxies are consistent with hosting a massive black hole at their formed intermediate-mass black holes in the early universe 

(M adau & Rees 2001). 

Numerical si mulations have demonstrated 

( Portegies Zwart et"al] l2004t iGiirkan et al.l l2004l iFreitag et"al] If the velocity dispersion - black-hole mass scaling rela- 



I2006h that intermediate-mass black holes can form in dense tion observed for supermassive black hole in galaxies (e.g. 
'Fen-arese & MerrittI l2000i iGebhardt et alll2000i iGiHtekin et al.l 

* Based on observations collected at the European Organization for 2009) holds at the lower end, intermediate-mass black holes 

Astronomical Research in the Southern Hemisphere, Chile (083.D- would be expected in systems with velocity dispersions between 

0444). 10-20 km/s like for globular clusters. Radio and X-ray detection 
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of gas in the central regions is also employed e ither to provide a 
black hole mass estimate or an upper limit (e.g. Maccarone et al.' 
,2005; Ulvestad et al. 2007; Bash et al. 2008; Cseh et al. 201ft 
Strader et al.l2012ft . While the current flux limits of lStrader et al.l 
(t2012) are impressively low, in order to provide an upper limit 
on a possible black hole mass, they must make various assump- 
tions about the gas accretion process. Some of the more uncer- 
tain assumptions include 1) the distribution of the gas, since a 
clumpy distribution will lead to time variability which signifi- 
cantly lowers the detection probability, 2) the efficiency of the 
accretion process, which may be particularly important if these 
systems are dominated by advection or convection, and 3) uncer- 
tainties in translating X-ray fluxes to bolometric fluxes to black 
hole masses. Thus, non-detection in X-ray and radio can be dif- 
ficult to interpret in terms of a black hole mass upper limit. 

So far the best candidates for hosting an IMBH are the most 
massive globular clusters in th e local group . One of the m is 
u) Centauri (NGC 5139) where 'Noyola et aL| (|2008l I20T 0I) and 
[jalali et al. (2012) detect the kinematic signafiire of a 40 000 Mq 
black hole based on radial velocities fro m integrated hght. The 
results were however challenged by ,van der Marel & AndersonI 
(1201 Ol) . who fin d a lower value of the black-hole mass via HST 
proper motions. IjaiaiT et al.l (l2012 f) perform N-body simulations 
and show that the current kinematic and light profile with re- 
spect to a kinematic center found in Noyola et al. (2010) are 
consistent with presence of a 5 x lO** Mq IMBH assuming a 
spherical isotropic model for this cluster. Another good candi- 
date is Gl in M31. It is the most massive globular cluster in the 
Local Group and is found to host a black hole of 20 000 
(eg. Gebhardt et alj 2005h by kin e matic measu r ements. In ad- 
ditioniPoolev & Rappaport(l2006l) : lKong(l2007h : IUlvestad et aP 
detected X-ray and radio emission at its center which is 
consistent with a black hole of the same mass. 

A third good candidate is NGC 638 8, a m assive globular 
cluster in our Galaxy. Lutzgendorf et'al] (|201 1) detected a rise 
in its central velocity-dispersion profile which is consistent with 
a black hole of ~ 20 000 Mq at its center The understanding of 
the formation and evolution of intermediate-mass black holes is 
crucial for the understanding of the evolution and formation of 
supermassive black holes. Seed black holes are needed in order 
to explain the fast formation process of these massive black holes 
which ar e obs e rved a t very high redshift, i.e. at early times in the 
Universe iFanI (120061) . Intermediate-mass black holes formed in 
globular c lusters an d accreted by their host galaxy could be these 
seeds (ee. lEbisuzaki et al. .2001: Tanaka & Haiman 2009). 

We chose to observe the globular cluster NGC 2808 as it 
shows a variety of interesting features. As known for some glob- 
ular clusters by now, NGC 2808 has multiple stellar populations. 
So far only one other cluster (a; Centauri) shows m ore than two 
distinguishable populations (eq. i Bedin et alJlIoo l. As observed 
by Piotto et al. (2007), NGC 2808 shows a triple main sequence, 
which indicates the existence of three sub-populations, all with 
an age of ~ 12.5 Gyr, but with different metallici ties. Also, its 
complex extended horizontal branch morphology (lHarrislll974l; 
iFerraro et al.lll990l) shows puzzling discontinuities in the stel- 
lar distribution along its blue tail. Maccarone & Servillat (2008) 
analyze deep radio observations of NGC 2808 and found no 
sources detected within the core radius. This places an up- 
per limit on a possible intermediate-mass black hole of 370 - 
2IOOM0 assuming a uniform gas density and the accretion rate 
to be different fractions of the Bondi rate. This limit can increase 
if one assumes even lower Bondi accretion rates or non-uniform 
gas content in the cluster 



Table 1. Propert ies of the globular c l uster N GC 2808 from the 
references: NG='Novola & Gebhardt (l2006l) . H= iHarrisI (Il996h 
and PM=Prvor & Mevlan ( 1993). 



I^alcllllC LCI 






RA (J2000) 


09h 12m 03s 


NG 


DEC (J2000) 


-64° 51' 49" 


NG 


Galactic Longitude I 


282.19 


H 


Galactic Latitude b 


-11.25 


H 


Distance from the Sun i?suN 


9.6 kpc 


H 


Core Radius rc 


12.4" 


NG 


Central Concentration c 


1.77 


H 


Heliocentric Radial Velocity V, 


101.6 + 0.7 km/s 


H 


Central Velocity Dispersion cr 


13.4 km/s 


PM 


Metallicity [Fe/H] 


-1.15 dex 


H 


Integrated Spectral Type 


F7 


H 


Reddening E(B-V) 


0.22 


H 


Absolute Visual Magnitude My, 


-9.39 mag 


H 



iNovola & Baumgardd (|201 ll) analyze N-body simulations of 
star clusters with and without central black holes. These reveal 
that the presence of an IMBH induces a shallow central cusp 
in the radial density profile. Hence, clusters showing shallow 
cusps are the best candidates for harboring an IMBH. Further, 
Miocchi (2007) investigate the effect of an IMBH on horizontal 
branch morphologies. A central black hole that strips enough 
stars of their outer envelope during close passages, could be 
one avenue for producing an extended horizontal branch (EHB). 
NGC 2808 displays a shallow cusp as well as an EHB, mak- 
ing it an excellent candidate for harboring a central black hole. 
Given its measured central velocity dispersion of cr = 13.4 km/s 
(iPrvor & MevlarLl993i) an extrapolation of the Mbh - cr relation 
predicts a black-hole mass of ~ 3 x IO^Mq. This translates into 
a radius of influence of 1 " - 2" at a distance of ~ 9.6 kpc, which 
would induce clear kinematic signatures inside the ~ 12" core 
radius. Further main characteristics of NGC 2808 are listed in 
table [U 

This work aims at investigating whether the globular cluster 
NGC 2808 hosts an intermediate-mass black hole at its center. 
We first study the light distribution of the cluster Photometric 
analysis, including the determination of the cluster center and 
the measurement of a surface brightness profile, is described 
in section |2] De-projecting this profile gives an estimate of the 
gravitational potential produced by the visible mass. The next 
step is to study the dynamics of the cluster Section |6] summa- 
rizes our FLAMES observations and data reduction and section 
|4] describes the analysis of the spectroscopic data. With the re- 
sulting velocity-dispersion profile, it is possible to estimate the 
dynamical mass of the cluster We then compare the data to Jeans 
models in section |5] Finally, we summarize our results, list our 
conclusions and give an outlook for further studies in section |6] 

2. Photometry 

The photometric data were taken from the archive of the Hubble 
Space Telescope (HST). They were obtained with the Wide Field 
and Planetary Camera 2 (WFPC2) in May 1998 (GO-6804, PI: 
F. Fusi Pecci) and are composed of a set of two exposures. The 
deep exposure dataset contains three images each in the filters I 
(F814W) and V (F555W) with exposure times of 120 and 100 s, 
respectively. In addition we use a set of two shallow images per 
filter with exposure times of 3 s in the I filter and 7 s in the V filter 
obtained in the same run. For both datasets the images cover 
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Fig. 1. The color-magnitude diagram of NGC 2808. Overplotted 
are the brightest stars identified in the ARGUS field of view (red 
circles), and the used template stars (green stars). 



the cluster center with the central 20" on the PC chip. The data 
were calibrated using the WFPC2-specific calibration algorithm, 
as retrieved from the European HST-Archive (ST-ECF, Space 
Telescope European Coordinating FacilitjQ). 

2.1. Color magnitude diagram (CMD) of NGC 2808 

The CMD is obtained using the programs daophotll, allstar and 
allframe by P. Stetson, applied to the HST images. F or a detailed 
documentation of these routines, see IStetsonl(ll987l) . These pro- 
grams are especially developed for photometry in crowded fields 
and are therefore ideally suited for the analysis of globular clus- 
ters. For the de scription of the individu al steps we refer to our 
previous paper (iLutz^endorf et al.ll201 lb . 

We obtain two catalogs with V and I magnitudes from the 
two datasets with different exposure times. We match and com- 
bine the catalogs using the routines CataXcorr and CataComb 
in order to obtain a complete star catalog over a wide magnitude 
range. All coordinates are transformed to the reference frame 
of the first I band image of the shallow exposure {u4fp010br). 
Figure [1] shows the final CMD with the positions of the bright- 
est stars in the ARGUS pointing and the spectroscopic template 
stars (see section |6]l overplotted. 



' Based on observations made with the NASA/ESA Hubble 
Space Telescope, and obtained from the Hubble Legacy Archive, 
which is a collaboration between the Space Telescope Science 
Institute (STScI/NASA), the Space Telescope European Coordinating 
Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre 
(CADC/NRC/CSA). 



2.2. Cluster center determination 

A crucial step is the determination of the cluster's center Precise 
knowledge of the center is important since the shape of the sur- 
face brightness and the angular averaged line-of-sight velocity 
distribution (LOSVD) profiles depend on the position of that 
center. Using the wrong center typically produces a shallower 
inner profile. Previous attempts have shown that the larger the 
core radius of the cluster, the more difficult it is to determine the 
exact position of the center For example, the exact position of 
the center of the globular cluster oj Centauri is still under debate 
and differs by up to 12" in different analyses (e.g. Novo la et al.l 
2010; Anderson & van der Marel 2010). In contrast, NGC 6388, 
with a core radius 10 times smaller than a> Centauri, turned out 
to have a well defined center when applying various method s to 
determine the center of the cluster (iLtitzgendorf et al.]|201 Ih . 

iNovo la & Gebhardt (l2006i hereafter NG06) determined the 
center of NGC 2808 to be ata = 09:12:03.09, 5 = -64:51:48.96 
(72000), with an uncertainty of 0.5", by minimizing the stan- 
dard deviation of star counts in eight segments of a circle. NGC 
2808 ha s a small core radius (12") in comparison to o) Centaur i 
(NG06; INovola et al.1 120101: lAnderson & van der Marell l20I0l) . 
In order to get an estimate of how accurately the center can be 
derived, we apply various routines to our catalog. 

The field of view of our dataset is very small because we are 
limited to the 34" x 34" field of the PC ship. In such a small 
field of view, it is difficult to determine the center given the large 
errors arising from the Poisson statistics, i.e. shot noise. Also, 
if the core is extended, we might not be able to see the stellar 
concentration decreasing. However, we compensate for this by 
using different techniques and by estimating the error from their 
scatter All techniques are applied for stars brighter than my = 20 
in order to account for the incompleteness effect for faint stars. 

The first technique uses isodensity contours as described in 
Anderson & van der Marel (2010). The field of view is divided 
into boxes of equal size of 100 x 100 pixels (4.6" x 4.6"). This 
size is the best compromise between having too few stars in each 
box and therefore large shot noise errors (boxes too small), and 
having not enough points and therefore a very noisy contour plot 
(boxes too large). The boxes contain about 100 stars on average. 
In each box, the stars are counted and the density is derived. The 
innermost isodensity contours, which were not disturbed by ge- 
ometrical incompleteness are fitted by ellipses and their central 
points are determined. From the average of these points and their 
scatter we determine the central position and its error. 

The second method is ve ry similar to the one desc ribed in 
McLaug hlin et alj ^2006) and iLiitzgendorf et all (1201 ll) . It uses 
a symmetry argument to determine the cluster center In our field 
of view a grid of trial centers is created, using a grid spacing of 
2 PC pixels (~ 0.1"). Around each trial center, a circle is traced 
and divided into wedges as shown in Figure [2] In order to in- 
crease the number of stars, we vary the size of the circle depend- 
ing on its position on the image. That is, the circle is always as 
big as the image borders allow it. The stars in each wedge are 
counted and the numbers compared to the opposite wedge. The 
differences in the total number of stars between two opposing 
wedges are summed for all wedge pairs and divided by the area 
of the circle. The coordinates that minimize the differences de- 
fine the center of the cluster We use different numbers of wedges 
from 4 to 16. We repeat the procedure by rotating the wedges 
so that their bisector is matching the x-y axis. This method is 
refined by comparing the cumulative stellar distribution of the 
stars in the opposing wedges instead of the star counts alone. 
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ISODENSITY CONTOURS PIE WEDGES LIGHT CONTOURS 

Density (arcsec"^) Summed Difference Mean Counts/Pixel 

2.3 3.6 4.9 6.2 7.5 8.8 10.1 431 369 306 244 181 119 56 1.1 1.5 1.9 2.3 2.7 3.2 3.6 




-10 -5 5 10 -4 -2 2 4 -10 -5 5 10 



arcsec arcsec arcsec 

Fig. 2. The method to determine the globular cluster center Left, ellipses fitted to isodensity contours. Middle panel, the contours of 
the cumulative pie wedges method. The wedges are shown in the upper right of the plot. The right panel shows the contours of the 
pixel analysis of the mapped image. Contours in blue mark areas of either high stellar density in the left panel, high symmetry (i.e. 
low differences between the wedges) in the central panel or high light density in the right panel. In every plot the triangle marks the 
center adopted by NG06 and the white star the newly derived center. 



The third and last method is based on the HST/WFPC2 im- 
age instead of the star catalog. This method computes the light 
distribution of the cluster For this we map out the brightest stars 
with circular masks in order to prevent a bias towards these stars. 
We divide the resulting image into boxes of 100 x 100 pixel and 
derive the mean value of the counts per pixel in each box. From 
this we compute a contour plot and fit ellipses to determine the 
center Despite our attempt to not get biased by the brightest 
stars, we find that in comparison to the previous two methods, 
the center obtained with the light distribution is shifted towards 
a clump of bright stars north-east of the center However, due to 
its large errors it does not bias the final center position. 

Figure|2]shows the contour plots of the different methods and 
Figure |3]presents a finding chart of our final center position. As 
a final result we obtain 

(xc,yc) = (344.4, 441 .6) + (9.9, 3.7) pixel (1) 

a = 09 : 12 : 03.107, Aa = 0.5" (J2000) (2) 
(5 = -64 : 51 : 48.45, A6 = 0.2" (3) 

which uses as reference image u4fp010br. We compared our 
result to the center obtained by NG06. The two centers are 
0.32" apart and thus coincide within the error bars of 0.5" 
(as determined by NG06 performing artificial image tests). 
This center also coinci des within 0.3" with the one derived by 
iGoldsburv et"an (1201 Ol) using ellipse fitting applied to the den- 
sity distribution on ACS/WFC images. 

2.3. Surface brightness profile 

The surface brightness (SB) profile is required as an input for 
the Jeans models described in the following section. As in 
iLiitzgendorf et al.l (1201 ll) we use a simple method of star counts 
in combination with an integrated light measurement from the 
WFPC2 image to derive the SB profile. For our analysis we do 
not use stars brighter than my = 17 to avoid contamination by 
very bright stars in the center The fluxes of all stars brighter 




Fig. 3. Finding chart for the center of NGC 2808. The magenta 
circle marks the center of NG06, the red ellipse the results from 
the pie wedges method, the orange ellipse is for the isodensity 
contour method and the yellow one shows the position and error 
of the method measuring the light distribution. The final adopted 
center and its error bar are displayed in green. 



than my - 20 are summed in radial bins around the center and 
divided by the area of the bin. In addition, the integrated light 
for stars fainter than my - 20 is measured directly from the HST 
image. Using the same radial bins as in the star count method, 
we measure the statistical distribution of counts per pixel exclud- 
ing regions containing stars with my < 20 by mapping out these 
stars in the image using a circular mask with a radius of 0.3" (6 
HST pixels). We use a robust bi-weight estimator to derive the 
mean counts-per-pixel for every bin. Finally, the flux per pixel 
is transformed into magnitudes per square arcseconds and added 
to the star count profile. Due to the small field of view of the 
PC chip we are only measuring points inside a radius of 10". 
The errors of our profile are obtained by Poisson statistics of the 
number of stars in each bin. With a linear fit to the innermost 
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log r (arcsec) 

Fig. 4. The surface brightness profile of NGC 2808. The red and 
the blue circles mark the measurements from the 2MASS im- 
age along the major and minor axis, respectively, as well as their 
MGE parametrization (dashed lines). The profile obtained from 
the HST star cat alog is shown in pu rple. Oveiplotted is the pro- 
file obtained bv lTrager et al.l (Il995b with a solid black line. 



points (r < 10") we derive a slope of the surface luminosity 
density I(r) o: r" of a - -0.16 + 0.08. This value is steeper (but 
consistent within the eiTors) than the slope of o' = -0.06 + 0.07 
derived by NG06. The final inner profile is listed in Table |2]and 
shown in Figure |4] 

For the outer regions we use the profile obtained by 
iTrager et al.l (1199 5'*) for our spherical models and a two- 
dimensional profile obtained from ground based (2MASS) im- 
ages (see Section I5.2l i for the axisymmetric models. The two 
2MASS images in J-band which cover the entire cluster were 
received from the public archive and combined to a single point- 
ing. Using routines provided in the package of the anisotropic 
Jeans models we obtain a two-dimensional surface brightness 
profile by fitting isophotes to the image. The final profile is con- 
structed by combining the inner profile obtained with the HST 
with the outer profile from the 2MASS image. Since the two im- 
ages are taken in different bands, it is necessary to scale the im- 
ages to a common flux. For this purpose we simply scale the data 
points of the 2MASS image to the HST profile. Here we neglect 
stellar population effects and assume a constant color within the 
cluster. Figure |4] shows the final combined profile with the mi- 
nor and major axis of the two-dimensional profile colored in blue 
and red, respectively. 



3. Spectroscopy 

The spectroscopic data were obtained with the GIRAFFE spec- 
trograph of the FLAMES (Fiber Large Array Multi Element 
Spectrograph) instrument at the Very Large Telescope (VLT) us- 
ing the ARGUS mode (Large Integral Field Unit). The obser- 
vations were performed during two nights (2010-05-05/06). The 
ARGUS unit was set to the 1 : 1.67 magnification scale (pixel 
size: 0.52", 14 x 22 pixel aiTay) and pointed to six different posi- 
tions, each of them containing three exposures of 600s with 0.5" 
dithering to cover the entire core radius. The position angle of 
the integral field unit remained at degrees (long axis parallel to 
the north-south axis) during the entire observation. 



Table 2. The derived surface brightness profile in the V-band. 
AVi, and AVi are the high and low values of the eiTors, respec- 
tively 



log r 


V 


A V 

LAV j 


A V;, 


[arcsec] 


[mag/arc sec^] 


[mag/arcsec^] 


[mag/arcsec-] 


-0.17 


15.00 


0.30 


0.24 


0.14 


15.07 


0.18 


0.16 


0.36 


15.27 


0.12 


0.11 


0.50 


15.18 


0.10 


0.09 


0.66 


15.32 


0.07 


0.06 


0.83 


15.39 


0.04 


0.04 


0.96 


15.48 


0.04 


0.04 


1.06 


15.58 


0.03 


0.03 


1.14 


15.75 


0.03 


0.03 



The kinematics are obtained from the analysis of the 
Calcium Triplet (~ 850 nm), which is a strong feature in the 
spectra. The expected velocity dispersions lie in the range 5-20 
kms"' and had to be measured with an accuracy of 1-2 kms"' . 
This implied using a spectral resolution around 10 000, avail- 
able in the low spectral resolution mode set-up LR8 (820 - 
940 nm, R = 10 400). 

We reduce the spectroscopic data with the GIRAFFE 
pipeline programmed by the European Southern Observatory 
(ESO). This pipeline consists of five recipes {gimasterbias, gi- 
masterdark, gimaster flat, giwavecalibra tion, giscience) which 
are described in Liitzgendorf et al.l(l201 ll). From the input obser- 
vations, the final routine giscience produces a reduced science 
frame as well as the extracted and rebinned spectra frame. At 
the end, the recipe also produces a reconstructed image of the 
respective field of view of the ARGUS observations. 

For sky subtractio n, we use the program developed by Mike 
Irwin and described in lBattaglia et al ] (l2008h . The program com- 
bines all 14 sky fibers using a 3-sigma clipping algorithm and 
computes an average sky spectrum. It splits the continuum and 
the line components for both the averaged sky spectrum and 
the object spectrum, using a combination of median and boxcar. 
The sky-line mask and the line-only object spectra are compared 
finding the optimum scale factor for the sky spectrum and the 
sky-lines are subtracted from the object spectra. The continuum 
is added back to the object spectra after subtracting the sky con- 
tinuum by the same scaling factor. 

As a next step w e use the program LA-Cosmic developed by 
Ivan Dokk um ( 200 1 ) to remove the cosmic rays from our spectra. 
In order to avoid bright stars dominating the averaged spectra 
when they are combined, we apply a normalization to the spectra 
by fitting a spline to the continuum and divide the spectra by it. 

4. Kinematics 

In this section we describe how we compute the velocity map in 
order to check for peculiar kinematic signals. Further, we mea- 
sure the velocity-dispersion profile which is used to fit analytic 
models, described in the next section. 

4.1. Velocity map 

To construct the velocity map, we use the relative shifts of the 
pointings to stitch them together and create a catalog with each 
spectrum correlated to one position in the field of view. The re- 
sulting catalog of spectra and their coordinates allow us to com- 
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bine spectra in different bins. The combined pointing contains 
54 X 44 spaxels and is cross shaped. We also reconstruct the 
ARGUS pointing on top of the HST image in order to obtain a 
direct connection between spaxel positions and our star catalog. 

Using the pena lized pixel-fitting (pPXF) program developed 
bv lCappellari & E msellem (2004), a velocity is derived for each 
spaxel. In regions of overlap the spectra are first averaged be- 
fore being analyzed. Figure |5] shows a) the field of view of the 
six ARGUS pointings on the HST image, b) the reconstructed 
and combined ARGUS image and c) the coiTesponding velocity 
map. For every ARGUS pointing, an individual velocity map is 
derived. The velocities of overlapping spaxels are in agreement 
with each other to 1 km s ' , which is expected as discussed in 
iLiitzgendorf et al.1 (1201 Ih . 

In order to test individual spaxels for shot noise, we apply the 
same routines as described in Liitzgendorf et al. (201 1 ). At every 
position of a star in the catalog, a two dimensional Gaussian is 
modeled with a standard deviation set to the seeing of the ground 
based observations (FWHM = 0.8") and scaled to the total flux 
of the star. We measure the absolute amount and fraction of light 
that each star contributes to the suiTounding spaxels. For each 
spaxel we then have the following information: a) how many 
stars contribute to the light of that spaxel, and b) which fraction 
of the total light is contributed by each star. The test shows that 
NGC 2808 is less concentrated in the center than NGC 6388 (in 
fact NGC 6388 has 1.5 times the central steUar density of NGC 
2808) and therefore fewer stars contribute to individual spaxel. 
Also, more spaxels are dominated by a single star by 60 % or 
higher. As in the case of NGC 6388, we map out spaxels in which 
either less than 10 stars contribute to the light or a single star 
contributes more than 60 % of the light. This leaves us with 1080 
spectra out of the 1514 spectra sample. 

The velocity template used in the kinematic analysis is im- 
portant to consider carefully given the strong changes in the in- 
trinsic line widths in a globular cluster. In order to find an op- 
timal velocity template we first collect all stars which dominate 
a spaxel by more than 80% via the shot noise routine described 
above. We plot these stars on our CMD in order to check for non- 
cluster members. We choose the five faintest stars and combine 
them after shifting them to the same velocity. We also try kine- 
matic fits with individual stars with high signal-to-noise from the 
upper giant branch. We do not find as good a fit to the integrated 
light using the stars from the upper giant branch as templates; 
this result is expected since the integrated light comes primarily 
from the fainter stars. The kinematics, however, are similar using 
both templates. 

The positions of the template stars are marked in the HST 
image (Figure 5, left panel) with magenta circles and in the CMD 
(Figure 1) with green stars. We also identify the brightest stars 
from the pointing in the CMD to make sure that none of the 
dominating stars is a foreground star (see Figure 1, red dots). In 
order to derive an absolute velocity scale, the line shifts of the 
template are measured by fitting a Gaussian to each Ca-triplet 
line of the template spectrum, and deriving the centroid. This is 
compared with the values of the Calcium Triplet in a rest frame 
and the average shift is calculated. The derived radial velocity is 
transformed to the heliocentric reference frame. This results in a 
template velocity of Vtemp = (122.3 + 2.3) kms"' . 

As a conspicuous feature in the velocity map we recog- 
nize two blue spots in the upper right, indicating two stars with 
high approaching velocities. Considering the velocity scale, also 
plotted in Figure |5] these features refer to velocities of -40 
and -45 kms ' relative to the cluster, respectively. This corre- 
sponds to 3.1 and 3.5 times the velocity dispersion (if one as- 
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Fig. 6. The combined spectra for the bins 1 , 3 and 6 from which 
the kinematic measurements are taken. 



sumes the value of iPrvor & Mevlanlll993h . High velocity stars 
have been discovered in only a few globu lar clusters up to now 
(iGunn & Griffinll979tlMevlan et alJl99ll) and require a detailed 
analysis in terms of membership, the tails of the velocity dis- 
tribution, and ejection mechanisms. We discuss the two high- 
velocity stars in a separate paper (Liitzgendorf et al., 2012 in 
preparation). 

4.2. Inner velocity-dispersion profile 

For the radial velocity-dispersion profile of NGC 2808, we bin 
the spectra in the following way. The pointing is divided into 
five independent angular bins, with radii of 2, 5, 10, 20 and 
28 ARGUS spaxels corresponding to 1.0", 2.6", 5.2", 10.4" and 
15.6" (~ 0.05,0.12,0.24,0.48,0.68 pc). We try diff^erent com- 
binations of bins and bin distances as well as overlapping bins 
and find no change in the global shape of the profile. In each 
bin, all spectra of all exposures are combined with a sigma clip- 
ping algorithm to remove any remaining cosmic rays. Velocity 
and velocity-dispersion profiles are computed by applying pPXF 
to the binned spectra using the same template as for the veloc- 
ity map. We compare the results of pPXF with a non-parametric 
fit. We obtain the best agreement of these two methods by fit- 
ting four moments of the Gauss-Hermite parametrization. The 
final parameters (V, cr, h^, h4) are transformed to the "true" mo- 
ments of the LOSVD (V, o", ^3, ^4) by applying equations (17) of 
van der Marel & Franx (1993). 

We estimate the radial velocity of the cluster in a heliocen- 
tric reference frame and the effective velocity dispersion cr^. 
We combine all spectra in the pointing and measure the veloc- 
ity relative to the velocity of the template. This value is cor- 
rected for the motion of the template and the heliocentric ve- 
locity and results in a value of Vr - (104.3 + 2.3) km s"' which 
agrees within the errors with the value from |Harris'('l996') = 
(101.6 + 0.7) kms ' . Th e effective velocity dispersion is de- 
rived using equation (1) in 'Liitzg endorf et al] (1201 II) and results 
in (Te = (13.4 + 0.2) kms ' . This is in very good agreement 
with the central velocit y dispersion cTc = (13.4 + 2.6) km s ' of 
iPrvor&MevTaiil (11993b . 

For the error estimation we run Monte Carlo simulations for 
each bin. From the routine described in section 14.1] we know 
how many stars contribute to what amount to each spaxel, and 
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Fig. 5. The velocity map of NGC 2808. Shown are the ARGUS field of view reconstructed on the WFPC2 image (left panel), the 
combined reconstructed ARGUS pointings (middle panel), and the resulting velocity map. Magenta circles mark the template stars 
used to derive the kinematics and the blue cross marks the center of the cluster. The white crosses on the velocity map mark the 
spaxel which are not used in deriving the velocity-dispersion profile. 



Table 3. The kinematics of NGC 2808 obtained from the VLT/FLAMES data. 
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how many spaxels are added up in each bin. Each of the stars 
in one bin is assigned a velocity chosen from a Gaussian ve- 
locity distribution with a fixed dispersion of 10 kms"' . Using 
our template spectrum we shift the spectra by their velocity and 
weight them according to their contribution before combining 
them in to one spaxel. The resulting spaxels are normalized, 
combined and the kinematics measured with pPXF (as for the 
original data). After 1000 realizations for each bin, we obtain 
the shot noise errors from the spread of the measured velocity 
dispersions. The eiTors for the velocity are derived by applying 
Monte Carlo simulations to the spectram itself. This is done by 
repeating the measurement for 100 different realizations, adding 
noise to the original spectra (see ,Cappellari & Emsellemii2004i, 
section 3.4). 

The resulting profile is displayed in Figure Q The innermost 
point drops down to a dispersion of 1 1 .4 km s ' , which is lower 
than the outermost point of the IFU data, but it is also severely 
affected by shot noise, as seen by its large error bar. In Table |3] 
we record the results of the kinematic measurements. The first 
column lists the radii of the bins. The following columns show 
the central velocities of each bin in the reference frame of the 
cluster, the corrected velocity dispersion & (Vrms), as well as 



the parameters from the Gauss-Hermite parameterization V, cr, 
h3 and h4. 



4.3. Outer kinematics 

In addition to the inner kinematics we use the dataset of 
Gebhardt et al. (2012, in preparation) for larger radii. This 
data was obtained in three epochs (1995, 1997, 1998) with the 
Rutgers Fabry Perot on the Blanco 4-m telescope at Cerro Tololo 
Inter-American Observatory (CTIO). The data set contains over 
3600 velocities of individual stars out to 4'. For reduction and 
more detailed information we refer to Gebhardt et al. (2012, in 
preparation). 

The Fabry-Perot velocities come from a v ery similar set 
of observations and reductions as presented in iGebhardt et al.l 
( 1997). These spectra are centered on a small region around the 
H-alpha absorption line, with absolute velocity calibration de- 
rived from comparison with published radial velocities. We ex- 
clude the Fabry-Perot data in the central regions where crowding 
is important; whenever measuring radial velocities of individual 
stars in clusters, one must be careful to limit contamination by 
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Fig. 7. Isotropic spherical Jeans models compared to the kinematic data of NGC 2808. Stars correspond to the inner kinematic 
profile obtained with integral-field spectroscopy and bullets represent the data obtained with the Fabry-Perot instrument. Shown are 
a model with a constant M/Ly over the entire radius of the cluster (black line) as well as the model where we fit an M/Ly profile 
(green line). The resulting M/Ly profile is shown in the right panel together with the constant value of the previous model. 



other cluster members, otherwise potentially biasing the velocity 
dispersion measurement. 

A closer look at the outer velocities shows a clear rotation 
at larger radii. We derive the rotation velocity and orientation 
in the following way: Different radial bins of 10", 30", 50", 70" 
and 90" are divided into 12 angular bins and the velocity of each 
of these wedges is derived. The rotation curve (the velocity as a 
function of position angle) is fitted by a sine function and orien- 
tation as well as rotation velocity is extracted for each of the bins. 
We obtain an average orientation angle of of maximum positive 
rotation of = 132° + 9° (measured from North to East). The ro- 
tation velocity increases from no rotation within the core radius 
to ~ 5 km s"' in the outskirts of the cluster. 

In order to obtain the velocity-dispersion profile, we apply 
the ma ximum likelihood method introduced bv iPrvor & Mevlani 
11991 . The iterative scheme used to solve the maximum like- 
lihood equations is very similar to that used by Gunn & Griffin 
( 1197 9) to fit the velocity scale parameter of King models. This 
approach is necessary especially when the uncertainties of the 
individual velocities are different. The data is divided into ra- 
dial bins and the velocity dispersion derived for each bin. Due 
to crowding in the cluster center, fainter stars get contaminated 
from background light which biases their measured velocities 
towards the cluster mean velocity. This results in a lower ve- 
locity dispersion for points measured in the central regions of 
the cluster. For this reason we only used measurements for radii 
larger than 50" (~ r/,). The resulting velocity-dispersion profile 
is shown in Figure |7] together with the points of the IFU mea- 
surements and the best fit Jeans model (see section|5]l. 



5. Dynamical Models 

After having extracted the velocity dispersion profile over a large 
radial range, the next step is to compare this data, together with 
the photometric profiles, to dynamical models. This section de- 
scribes the different types of models which we compare to our 
data. We start with a spherical isotropic Jeans model, increase 
the complexity of the models to axisymmetric models, and ra- 
dial varying M/Ly profiles. 



5. 1 . Isotropic spherical Jeans models 

Assuming spherical symmetry is a good approximation for most 
globular clusters and a valid first order assumption for NGC 
2808. To compute the models, we use the Jeans Anisotropic 
multi-Gaussian expansion (JAM) dynamical models implemen- 
tation for stellar kinematics of spherical and axisymmetric 
galaxies developed by Cappellari (2002, 2008) 0. The rou- 
tines take a one dimensional surface brightness profile as an 
input an d use t he mu lti-Gaussian expansion (MGE) method 
Emselle m et alj d 19941) to fit and deprojects the light profile. 
This is then fed into the spherical Jeans equations and a second- 
moment profile (Vrms = -^(cr^ + V^^,), hereafter referred to as 
velocity-dispersion profile) is computed. The modeled velocity- 
dispersion profile is scaled by a constant factor to fit the kine- 
matic data. This scaling factor is adopted as the global M/Ly 
value. In Figure |7] the black line in the left panel represents the 
spherical Jeans model with a M/Ly = 1.4. 

The comparison of the spherical model with our data in 
Figure |7] shows already good agreement. However, a constant 
M/Ly profile over a larger radial range is not a good assump- 
tion for a globular cluster. N-body simulations have shown that 
the M/Ly increases for larger radii due to mass segregation and 
migration of low-mass stars towards the outskirts of the clus- 
ter. The underestimation of the model for the outer points might 
be caused by this effect. Therefore, it is necessary to allow for 
a varying M/Ly profile in the models. For this, we apply two 
methods: the first one is to let the model fit the M/L profile 
to the data. This can be done by multiplying each of the MGE 
Gaussians of the fitted surface brightness profile with different 
factors p ( fWilliams et alJI2009l) . These factors are varied over 
a physical range (p e (0, 3}). For every combination, the Jeans 
model is computed and the quality of the fit calculated via a least 
square statistic. With this technique, we find the M/Ly profile 
which best reproduces the data. The disadvantage of this method 
is the degeneracy of the problem. There are many combinations 
of the factors p that return similar quality of the fit to the kine- 
matic data. In addition, without constraining the parameters of 
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the M/Ly profile to a physical limit, the best fit would favor un- 
realistic high values in the outer regions to fit the flat velocity- 
dispersion profile beyond 100". The resulting M/Ly profile is 
shown in the right panel of Figure |7] as a green line. 

The second method takes an already existing M/Ly profile 
obtained from N-body simulations and is explained in the next 
Section. 

5.2. Isotropic axisymmetric Jeans models 

The Fabry-Perot dataset shows a clear rotation in the outer re- 
gions with rotation velocities up to 5 km s ' . The shape of NGC 
2808 reveals the rotation character of that cluster with a flatten- 
ing of e = 1 -Z7/a = 0.12 (IWhite & Sh awl 1987). For that reason 
we compute axisymmetric Jeans models in addition to the spher- 
ical models described in the previous section. These models as- 
sume an axial symmetry of the velocity ellipsoids rather than a 
spherical symmetry and are ideally suited for rotating systems. 
By a pplying the general axisymmetric Jeans equations (3) and 
(4) in lCappellaril (l2008h . the Jeans models provide good descrip- 
tions of the two-dimensional shape of the velocity (V) and its 
velocity-dispersion (Vrms), once a surface brightness profile is 
given. 

Instead of using an one-dimensional surface brightness pro- 
file the axisymetric models (also present in the JAM package) 
perform a two-dimensional fit to the surface brightness obtained 
from an image of the object. For globular clusters, determining 
a 2D surface brightness profile is challenging since the individ- 
ual stars are resolved and the routine has difficulties in fitting 
isophotes to a discontinuous image. Images with a lower resolu- 
tion are better suited for this purpose. We use a combination of a 
ground based image in J-band taken from publicly available data 
of 2MASS and of our inner one-dimensional surface brightness 
profile obtained from the HST image (see Section l27Jt . Since the 
high resolution of the HST does not allow a two-dimensional fit 
on the image itself and because we do not measure any rota- 
tion within the core radius, we assume spherical symmetry in- 
side a radius of r ~ 15"(the size of the combined ARGUS field 
of v iew). This assumption is further supported by the results of 
iLiitz gendorf et al. (2011) where we found that any anisotropy in 
a globular cluster is smoothed out within a few relaxation times 
inside the core radius of the cluster 

For the outer region of the cluster the routine find, 
galaxy.pro, which is also included in the JAM package, mea- 
sures an ellipticity of e = 1 - b/a = 0.11 and an orientation of 
6 - 104°. This orientation differs from the one that we obtain 
from the kinematic studies (0 - 132°) and from the results of 
[white & Shawl ( 1987) of 9 ~ 121°. Because of the large uncer- 
tainties of the photometric measurement (due to shot noise of 
single stars), we decide to adopt the orientation from the kine- 
matics as our final result - 132° + 9°. The next step is to de- 
termine the two-dimensional surface brightness profile from the 
2MASS image. This is done by the routine sectors- photome- 
try.pro. The routine performs photometry of an image along sec- 
tors equally spaced in angle. The result is converted into magni- 
tudes per squarearcsecond and combined with the inner surface 
brightness profile obtained with the HST. 

The JAM code calculates a two dimensional Vrms rnap by 
using the surface brightness profile. To compare with our one 
dimensional kinematic profile we extract the profile along the 
major axis and the minor axis. In addition, we derive a third pro- 
file by averaging the two-dimensional Vrms map over concentric 
radial bins and call it the averaged profile. Figure |8]shows these 



three profiles on top of the data. The models differ only in the 
outer region since the inner part is assumed to be spherical. 

In order to use the most realistic M/Ly profile, we use an 
M/Ly profile obtained by running N-body simulations and com- 
paring them to our data points using the N-bo dy code NBODY6 
(Aarseth 1999). Starting from a Kin^ (Il962h model the simu- 
lations run with varying initial central concentrations c, half- 
mass radii r/,, and IMBH masses with = 50 000 parti- 
cles. We use a metallicity of [Fe/H] = -1.18 and model a 
stel lar evolu t ion ac cording to the stellar evolution routines of 
Hu riev et al.1 (l2000l). The simulations start with stars distributed 
according to a Kroupa (2001) mass function in the mass range 
0.1 < ih/Mq < 100. Binaries can form during the evolution of 
the cluster, but are not included primor dially. For more inf orma- 
tions about these simulations we refer to lMcNamara et all(l2012h 
and Baumgardt & Makino (2003). We use a grid of computed 
models with different initial conditions and find the model which 
fits the observed kinematic and light profile the best. The M/L 
of this model is then computed as a function of radius by only 
using the brightest stars. This allows an independent determina- 
tion of the M/L profile of NGC 2808. This is fed into the Jeans 
model, which parametrizes the input M/L profile with an MGE 
fit and applies it to the resulting velocity-dispersion profile. The 
lower right panel of Figure |8] shows the used M/Ly profile from 
the N-body simulations (black line) as well as its parametriza- 
tion with the MGE Gaussians (green line). Also shown is the 
resulting model compared to the data (lower left panel). Model 
provides a worse fit then the model with a constant M/Ly pro- 
file. This comes probably from overestimating the M/Ly profile 
in the outer regions of the cluster 

Both methods of deriving an M/Ly profile show similar re- 
sults. Figure |8]in the upper panel shows the Jeans model for the 
major axis, minor axis and circular averaged on top of our data 
points. The model of the major axis seems to reproduce the data 
best. This might be due to the fact that the data set of the Fabry- 
Perot observations is asymmetric and spatially biased towards 
the major axis. The right panel of the figure shows the depro- 
jected M/L profile from the N-body simulations (black line) to- 
gether with the fit of the M/L profile from the Jeans model (red 
line). The M/L profile in the lower panel of Figure |8] has an 
interesting shape. The steep rise at the center implies a high con- 
centration of stellar remnants and therefore an advanced stage of 
mass segregation. Beyond the half-light radius (~ 48") the M/L 
profile rises again, which can also be explained by the process 
of mass segregation. Low mass stars move towards the outer re- 
gions while the cluster evolves. That explains the higher M/L 
ratio in the outskirts of the cluster 

The next step in terms of modeling is to include a central 
black hole in our Jeans models and to test if we obtain a better fit 
to our data. Figure |9] shows the result of these models. For both 
cases, with and without a M/Ly profile, we compute models with 
black-hole masses between Mbh - and Mbh = 8x10"' Mq. 
The black solid line shows the best fit of the model with a zero 
mass black hole. Higher black-hole masses predict higher cen- 
tral Vrms than seen in our data. The curves in the right panel 
of the figure shows this result. The blue shaded areas define the 
Icr limit of the best fit. For the model with the constant M/Ly 
profile, no black hole is needed in order to reproduce the data 
with a Icr uncertainty of M, = 1 x IO^Mq. The second model, 
however, exhibits a flatter profile and therefore is better fit by 
a model with a non-zero black-hole mass. The best fit model 
returns a black-hole mass of M, = 4 x 10'* Mq and a Icr up- 
per limit of M. ~ 6 X IO^Mq. Also the total M/Ly values of 
both models differ The model with the constant M/Ly results 
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Fig. 8. Axisymtnetric Jeans models for NGC 2808 for a constant (upper panel) and a varying M/Ly profile (lower panel). The used 
M/Ly profiles/values are shown in the right panels. The black line marks the inserted profile from N-body simulations and the green 
line its parametrisations with the MGE Gaussians. 



in a total value of M/Ly - 1.6Mq/Lq while the model with the 
M/Ly increasing in the outskirts of the cluster predicts a M/Ly 
of ~ 2.IM0/L0. 



We run Monte Carlo simulations on both the surface- 
brightness profile and the velocity-dispersion profile in order 
to get an estimate of how much the individual errors from the 
two profiles influence the result. We find that by only chang- 
ing the inner surface-brightness profile, 90 % of the cases re- 
sult in the same black-hole mass as the fit with the original light 
profile. This indicates that the error on the black-hole mass is 
only dependent on the uncertainties of the velocity-dispersion 
profile. We therefore adopt the 95% confidence upper limit of 
M. < 1 X 1O'*M0 of the model with the varying M/Ly profile. We 
note that the model with the constant M/ Ly fits the data slightly 
better and results in a lower upper Umit of M. < 3 X IO^Mq. 
However, a constant M/Ly profile is very unrealistic for a glob- 
ular cluster like NGC 2808 and therefore the results from the 
model with the M/Ly profile more reliable. From the Monte 
Carlo simulations we also derive a Icr error on the total mass 
and the global mass to light ratio of ~ 10 %. 



6. Summary and Conclusions 

We examine the central kinematics of the globular cluster NGC 
2808 in order to constrain the mass of a possible intermediate- 
mass black hole at its center With a set of HST images, the 
photometric center of the cluster is determined. Furthermore, a 
color magnitude diagram as well as a surface brightness profile, 
built from a combination of star counts and integrated light, are 
produced. The spectra from the VLT ground-based integral-field 
unit ARGUS are reduced and analyzed in order to create a veloc- 
ity map and a velocity-dispersion profile. We derive a velocity- 
dispersion profile by summing all spectra into radial bins and ap- 
plying a penalized pixel fitting method. In addition, we use radial 
velocities measured with the Rutgers Fabry Perot by Gebhardt et 
al. (2012, in preparation) to derive the velocity-dispersion profile 
in the outer regions. 

We compare the data to spherical and axisymmetric isotropic 
Jeans models using different M/Ly profiles and black-hole 
masses. We include a radius dependent M/Ly profile obtained 
from N-body simulations in order to account for mass segrega- 
tion in the cluster This does not improve the fit of the velocity 
dispersion profile, but still coincides with the data within the er- 
ror bars. The best fitting Jeans model is the axisymmetric case 
with a constant M/Ly profile and no black hole. The 95% error 
of that fit predicts an upper limit of M, = 3 x lO^M© on the 
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Fig. 9. Axisymmetric Jeans models with different black-hole masses for NGC 2808. Upper panel: The model with the constant 
MjLv (same as in Figure[8ll. Lower panel: The model with an MjLy profile. The jchfi values as a function of black-hole mass are 
shown in the right panels. Blue shaded areas mark the black-hole mass range where l^x^ < thus the Icr error. 



black-hole mass. However, assuming a non constant M/Ly pro- 
file the model upper limit on the black-hole mass increases to 
M. = 1 X IQ'^Mq. 

Our result on the upper limit of the black-hole mass in 
N GC 2808 is higher with the re sults of the radio observations 
of|M 

accarone fe Servillad dloosl) and their upper limit of 370 
- 2100 Mq. With the uncertainties in gas content and accretion 
rates, however it is plausible for the limit derived by the radio ob- 
servations to increase. We stress that the result is dependent on 
the choice of the M/Ly profile and needs to be treated carefully 
since different M/Ly profiles bring different results on black- 
hole mass and total mass of the cluster. 

Our derived mass to light ratio of M/Ly - (2.1 + 
0.2) Mq/Lq is higher than the M /L of ~ 1.3 derived by 
iMcLaughhn & van der Marel (l2005h . This results from the fact 
that ou r total luminosity is lower than the one derived by Harris 
(11996 *) by ~ 22% which results from the slightly steeper drop 
of the s urface brightness pr ofile compared to the profile ob- 
tained bv lTrager et al](ll995l) in the outskirts of the cluster. Also 
the total mass of Mtot - (8.2 + 0.8) x 10 ^ Mp is higher than 
the total mass derived with t he va lues of iHarrisI (119961) and 
iMcLaughUn & van der Marel (l2005l) of Mtot ~ 6.4 x IO^Mq 
which results from the shape of our M/Ly profile. 



So far, Jeans models allow us a crude first guess on the dy- 
namic state of a globular cluster. Nevertheless, the result of the 
black hole-mass depends strongly on the M/Ly profile used. We 
find that using different profiles with lower M/L values at the 
center but still high values in the outskirts result in models which 
fit the data with a black-hole mass up to 5000 Mq. However, due 
to the high degeneracy we encounter when fitting the M/L pro- 
file with Jeans models, we have to choose the profile which is 
derived from N-body simulations. Therefore, it is crucial to run 
specific N-body simulations for all our globular cluster in order 
to get an accurate M/Ly profile and constraints on the anisotropy 
and mass segregation. This is an important factor especially for 
mass segregated clusters such as NGC 2808. 

The dynamical models presented here include Jeans 
isotropic modeling and comparison with N-body simulations. 
We do not present results from orbit-based models. These 
axisymmetric models are significantly more general than the 
isotropic models since they have no assumption about the veloc- 
ity anisotropy. They are also more general than the N-body sim- 
ulations, since the N-body models rely on a limited set of initial 
conditions; the orbit models encompass all available phase space 
configurations, at the expense of producing a dynamical model 
that does not take the evolutionary processes into account. Thus, 
the orbit models will provide larger uncertainties, and hence a 
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large upper limit on the black hole mass. For this analysis, how- 
ever, we rely on the N-body simulations which should be a fair 
representation of the current dynamical state of the cluster. 

The study of black holes in globular clusters has drawn the 
attention of the astronomy community. Not only radial veloci- 
ties are observed and analyzed, also, observations and analysis 
of proper motions for many clusters are in progress. A desired 
future project would be the combination of all these data sets 
and a detailed analysis via N-body and orbit based models which 
would allow a deeper insight into the dynamics of globular clus- 
ters, revealing their hidden secrets. 
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